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Abstract. We consider the free fall of slender rigid bodies in a viscous incompressible fluid. We 
show that the dimensional reduction (DR), performed by substituting the slender bodies with 
one-dimensional rigid objects, together with a hyperviscous regularization (HR) of the Navier- 
Stokes equation for the three-dimensional fluid lead to a well-posed fluid-structure interaction 
problem. In contrast to what can be achieved within a classical framework, the hyperviscous 
term permits a sound definition of the viscous force acting on the one-dimensional immersed 
body, and global-in-time existence and uniqueness of a solution can be proved. Those results 
show that the DR/HR procedure can be effectively employed for the mathematical modeling of 
the free fall problem in the slender-body limit. 

1. Introduction 

Composite systems where microscopic bodies move in a viscous fluid are ubiquitous in both 
biological and technological contexts. The shape of such immersed bodies often displays some 
slenderness: one can encounter flat bodies, elongated bodies, somehow entangled filaments, and 
often simply point-like particles. For both mechanical and computational reasons it is tempting to 
approximate slender three-dimensional bodies with lower-dimensional objects, and a slendcr-body 
theory, especially for low-Rcynolds-number flows, has been widely developed (see e.g. [1, 2, 13, 
14, 12, 11]). Nevertheless, such a theory often relies on asymptotic expansions and sometimes 
crude approximations, mainly because of the difficulty in dealing with complex geometries and 
of the lack of a sound mathematical model for the interaction between three-dimensional fluids 
and lower-dimensional immersed bodies. The key obstacle to the development of such models is 
discussed in Remark 1. 

Here we follow the spirit of a new slender-body theory proposed in [6], which, in contrast to 
previously developed alternatives, might be viewed as a modeling scheme, based on the concepts of 
dimensional reduction and hyperviscous regularization. The geometry of slender bodies immersed 
in a viscous fluid is approximated by lower-dimensional objects, while adding a hyperviscous term 
to the flow equation. The hyperviscosity, an additional parameter entering the equation, is given 
by the product of the ordinary viscosity with the square of a length that replaces the characteristic 
size of the body along the dimensions that vanish in the slender-body limit. 

The classical Navicr-Stokes equation for incompressible Newtonian fluids reads 

du . , . 

p-^ + p[u ■ V)u + \7p - pAu = pb , 

where p is the pressure field, u is the divergence-free velocity field, p > is the constant and 
homogeneous mass density, p > is the dynamic viscosity and pb is a volumetric force density. It 
is well-known that the regularity of solutions for the three-dimensional Navier-Stokes equation is 
still an open problem, and various modifications of that equation with better regularity theories 
have been analyzed. Among those, the hyperviscous regularization (see, for example, [15], Chap. I, 
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Remarque 6.11) consists of adding a term proportional to AAm to the equation. For this modified 
equation, 

p— + p(u ■ V)it + Vp - pAu + CAAm = pb , 
at 

where £ > is the hyperviscosity, existence and uniqueness of regular solutions have been estab- 
lished. 

Despite the mathematical appeal of the hyperviscous regularization, it is generally hard to assign 
a relevant physical meaning to £, other than that of a higher-order dissipation coefficient. In a 
series of papers [3, 16, 8, 5], different contributions to ( associated with dissipation functional are 
introduced and analyzed. Here, as stated above, we assign to C a geometric, rather than dynamical, 
meaning, by introducing the effective thickness L > of the lower-dimensional objects, and setting 
£ = pL 2 , so that the hyperviscous flow equation becomes 

p— + p(u-V)u + S7p- pA(u- L 2 Au) = pb. (1) 

It has been shown in [7] that the flow generated by one-dimensional rigid bodies, moving 
with a prescribed velocity in a fluid governed by the hyperviscous flow equation, can be uniquely 
determined. The purpose of the present note is to establish the existence and uniqueness of 
solutions for the free fall of a slender rigid body E in a hyperviscous fluid. In the problem at 
hand, the velocity of the body is no longer prescribed, and it has to be determined, together with 
the fluid velocity, by suitably coupling the flow equations with the equations describing the motion 
of the rigid body. The free fall of three-dimensional rigid bodies in a Newtonian fluid has been 
widely studied; an excellent review of this topic can be found in [4], and we will follow similar 
arguments in our analysis. The main novelty here is, of course, the gap between the dimensions 
of the body and those of the fluid, which can be handled by a suitable definition, presented in 
Section 3, of the hydrodynamic force exerted by the hyperviscous fluid on the rigid body. 

2. Formulation of the free fall problem 

The free fall problem is characterized by the fact that the rigid body is immersed and dropped 
from rest in an otherwise quiescent fluid filling all of space, and gravity is the only external force 
acting on the system. We represent a rigid body as a connected and bounded closed subset E 
of M 3 which is a finite union of images of [0, 1] through a C^-diffcomorphism. In this way, the 
one-dimensional Hausdorff measure of E, H 1 (E), is finite and positive. It is convenient to write 
the problem in a co-moving frame, with origin at the center of mass c(f) of E. Denoting by y the 
position of a point in the original incrtial frame, and by x its position in the co-moving frame, we 
know that, at any time t > 0, 

as = QT(t)(v-c(t)), 

where Q(t) is an orthogonal linear transformation for any t > 0, with Q(0) = 1. If the velocity of 
the center of mass and the spin of the rigid body in the inertial frame are denoted by 77(f) and 
J?(f), respectively, so that 

v(t) = 77(f) + n(t) x{y- c(f)) 

is the velocity, in that frame, of any point y belonging to the rigid body, then their expression in 
the co-moving frame is given by 

£(t) := Q T (f)77(f) and w(t) := Q T (t)f2(t) , 

respectively, and the rigid velocity field v is transformed into 

17(f) := £(*) +w(t) XX, 

where x denotes the coordinates of y in the co-moving frame. 

We also note that the gravitational acceleration vector g (constant in the inertial frame) is 
represented, in the co-moving frame, by G(f) := Q T (t)g, which is easily seen to satisfy the ordinary 
differential equation 

dG „ 

- = Gx.. (2) 
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As customary when studying flows past rigid bodies, the velocity field u that we consider is the 
so called disturbance field, which is the difference between the actual flow and the flow at infinity, 
both seen in the co- moving frame. Since the flow at infinity is — C/ (that is minus the extension to 
all of the fluid of the motion of the immersed object), the representation of the fluid flow in the 
co-moving frame is given by u — U. 

The continuity and flow equations for an incompressible (disturbance) velocity field u(x, t) and 
pressure field p(x,t) defined on (R 3 \ E) x [0, +00), become 

divw^O, (3) 



p + [(u - U) ■ V}u + u xu^j = divT(tt,p) + pG , (4) 

where T(u,p) denotes the Cauchy stress tensor. Note that, thanks to its frame indifference prop- 
erties, T retains the same functional dependence on the velocity field seen in both the inertial 
frame and in the co-moving one. 
We have also a decay condition 

lim u(x, t) = , (5) 

I as I — >oo 

and the adherence to the rigid body, given by 

u(x,t) = U(x,t) on E x [0, +00) . (6) 
The equations of motion for the rigid body in the co-moving frame are then 

+ muj x £ = m e G + f (7) 



and 



+ u> x (Ju>) = -m c (r xG)+t, (8) 



where m is the total mass of the rigid body, m e is its effective mass, m c is the complementary 
mass, J its inertia tensor, r is the position of the centroid*- 1 - 1 of E in the co-moving frame, and 
/ and t are the total hydrodynamic force and torque, respectively, exerted on the rigid body as 
a consequence of the fluid flow u and pressure p. The proper definition of / and t is discussed 
in Section 3. In view of the interpretation of one-dimensional bodies as representations for real 
three-dimensional objects, to properly account for Archimedean forces which would vanish in the 
slender-body limit, we have to consider from the very beginning m e = m — m c as the difference 
between the real mass of the object and the mass m c of a portion of fluid occupying the real 
volume of the object, and we must add the terms proportional to m e and m c in equations (7) 
and (8), respectively. 

The whole set of equations (2)-(8) represents the differential problem associated with the free 
fall of a rigid object E in an incompressible fluid. The physical properties of such a fluid are 
encoded in the Cauchy stress tensor T, and the peculiar form we are going to use throughout the 
present paper is discussed in Section 3. 

It is convenient to consider the non-dimensional form of problem (2)-(8); by choosing suitable 
reference speed W and length d, we can switch to non-dimensional quantities according to 

x tu u pd m 

d pd 1 W fxW pd A 



WThe centroid of a rigid body coincide with its center of mass when the body has a uniform mass density; in 
the latter case r = 0. 
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obtaining 



div« = 0, (9) 
du 

— + Re{[(u - U) ■ V]u + u xu} = divT{u,p) + G, (10) 

lim u(x,t) = 0, (11) 

\x\— >oo 

u(x,t) = U(x,t) on E x [0, +oo) , (12) 

dt 
ice 
dt 
IG 
dt 



+ Re(mu> x g) = m e G + f , (13) 
J^ + &[ux(Ju)] = -m c (rxG)+t, (14) 
^- = Re(Gxuj) 1 (15) 



with initial conditions 

u(x,0)=£(0) = a;(0) = 0, G(0) = g . (16) 

In the foregoing equations Re = pWdj ix is the Reynolds number and every quantity has to be 
understood as non-dimensional. Moreover, a suitable choice of the reference quantities ensures 
that |G(t)| = 1 for any t > 0, that is G takes values on the sphere S 2 . The low- Reynolds- number 
approximation of the differential problem is obtained by neglecting the terms proportional to Re 
in equations (10), (13), and (14), whereas equation (15) remains unchanged, since it represents a 
geometric constraint which holds for any non- vanishing value of Re. 

3. The viscous force acting on a slender body 

The constitutive theory for non-simple fluids leading to a hyperviscous flow equation has been 
developed in [3, 16, 8]. It offers a number of possible choices for the terms to be included in T, in 
addition to those of Newtonian fluids. Here we make a somewhat minimal assumption, obtaining 
a fluid which is quasi-Newtonian, while being able to adhere to lower-dimensional objects. As 
explained in [5], we can simply take 

T(u,p) := -pl+ (Vu + Vu T -£ 2 VAu) , (17) 

where we used the non-dimensional effective thickness £ = L/d. That choice for T has also the 
feature of introducing only one new parameter, £, to which we have already assigned a geometric 
meaning. It is straightforward to check that T, defined as in (17), enjoys the standard symmetry 
and frame indifference properties. 

Let C be any closed bounded subset of E 3 and r > 0. The set 

V r (C) := {x e R 3 : d(x, C) < r) , 

where d(x, C) denotes the distance of x from C, is the r-neighborhood of C. We define the total 
hydro dynamic force, due to the fluid velocity and pressure field (u,p), acting on the slender body 
E as 

f(u,p) := lim / T(u,p)n, (18) 

where n denotes the unit outer normal to <9K-(E), and T is the stress tensor defined in (17). Notice 
that, thanks to the regularity of E, there exists always f > small enough so that V r (T,) has a 
Lipschitz boundary for any r < f. 

Let us check that the limit in (18) is indeed well-defined. We consider a ball Br centered at the 
origin and with very large radius i?, which contains E. According to equation (1), the term divT 
is required to balance both the inertial term and the external forces acting on the system. In the 
free fall problem considered in this paper, the only external interaction is gravity and, due to the 
dimensional reduction we perform on the rigid body, it is represented by a measure whose singular 
part is concentrated on E. We can see this by noting that the mass density per unit volume must 
diverge on E to give a non-zero weight to a body with vanishing volume. 
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Hence we expect also div T to be a measure and, denoting by a (?) its absolutely continuous 
(singular) part with respect to the three-dimensional Lebesguc measure, we expect the support of 
? to be contained in E and that divT = a in R 3 \ E. Under these assumption we have, for any 
r > small enough, 



lim / div T = lim / a = / a = div T , 

JB n \V r (T,) JB R \V r (T.) JB n \S JB a \T, 

by applying Lebesgue's theorem, since a GL'fi 3 ); then, by the Diver gence theorem, 



dV r (T.) JdB R J B R \S 



/=lim / Tn= / Tn- divT, (19) 



where n is always the outer normal. Since the right-hand side of (19) is independent of r, the 
left-hand side is well defined. It is important to stress the fact that, if ? were absent, then the 
integral over Br \ E of div T would be equal to the integral over all of Br and / would simply 
vanish. 

In a similar fashion, we define the total hydrodynamic torque acting on E, due to the fluid 
velocity and pressure field (u,p), as 

t(u,p) := lim / xxl(u,p)n. (20) 



9V r (E) 



4. Steady free fall 

In what follows we prove the existence of a solution for the steady version of the differential 
problem introduced in sec. 2, which consists in finding (u,p) £ if 2 (R 3 ;R 3 ) x ff _1 (R 3 ;R) and 



C,uel 3 , and geS 2 such that 

divw = 0, (21) 

Re{[(u - U) ■ V]u + lo x u] = divT(w,p) +g in R 3 \ E , (22) 

u{x) = £ + w x x onE, (23) 

Re{muj X £) = m e g + / , (24) 

Re[u) x (Jw)] = — m c (r xG)+t, (25) 

9xu = 0. (26) 



We will first treat the low- Reynolds-number version of problem (21)-(26), thus reducing the non- 
linearity to the sole term gxo). The linear nature of the flow equations in the low- Reynolds- number 
limit allows for a finer analysis of the relation between the shape of the objects and their steady 
motion, which is developed in [9]. Here we only present an existence result, studying then the full 
nonlinear problem, which is the main topic of the present paper. 

4.1. The low-Reynolds-number limit. The differential problem considered in this section con- 
sists in finding (u,p) £ # 2 (R 3 ;R 3 ) x fl r - 1 (R 3 ;M) and ? ,ue R 3 , and g <E S 2 such that 



div« = 0, (27) 

Vp + C u = g in R 3 \ E , (28) 

u(x) = £ + uj x x onE, (29) 

m e g = f , (30) 

m c (rxg)=t, (31) 

gxu) = 0, (32) 



where we set Cu := — Am + ^ 2 AAm. 

To prove the existence of a solution of the steady free fall problem, we will introduce some 
auxiliary problems, which are well-posed by virtue of the following result. 
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Lemma 4.1. Given £ G R 3 and uj G R 3 , t/iere exists a unique solution (h,p) of the problem 

div h = m R 3 , 

Vp + C h = m R 3 , (33) 
ft, = £ + u> x cc on S , 

belonging to the space if 2 (R 3 ;R 3 ) x ii _1 (R 3 ;R). Moreover, such a solution is of class C°° on 
R 3 \E. 

Proof. Consider the space 

H := {v G i? 2 (R 3 ;R 3 ) : divv = and v = £ + w X x on E} , 

which is closed and convex in if 2 (R 3 ; R 3 ), since the latter embeds in a space of Holder continuous 
functions. Then, the velocity field h can be found by minimizing on H the functional 

F{v):= \ f (2|SymVt;| 2 +i 2 |Au| 2 ) , 

where Sym Vd := (Vv + \7v J )/2. Being T a convex functional, h is unique, and the pressure 
field p can be recovered as the Lagrange multiplier of the divergence- free constraint. The same 
technique used to study the regularity of the solution for the classical Stokes equation can be 
applied here to prove that the solution is of class C°° on I 3 \ S. □ 

Remark 1. In view of the proof of Lemma 4.1 it is clear that, for h to truly satisfy the adherence 
condition on E, the latter must have positive i? 2 -capacity. That property is granted by the embed- 
ding of if 2 (R 3 ;R 3 ) into a space fo Holder continuous functions. On the contrary, when treating 
the classical Stokes equation, the natural functional space for the velocity field is 7? 1 (R 3 ;R 3 ). 
Now, the if ^capacity of S vanishes, allowing us to arbitrarily change the value of the velocity 
field on E, so that the imposition of an adherence condition does not really affect the flow. Indeed 
it is immediate to see that, given £ G R 3 and us G R 3 , the solution u G H 1 (M 3 ; R 3 ) of the problem 

div u = in R 3 , 

Vp - pAu = in R 3 , 

u = £ + uj x x on E , 

is such that 

M 2 = o, 



and hence u = almost everywhere in R 3 . 

Consider now the solutions (7i (i) ,f>W) an d (H®,P&) (i = 1,2,3) of 

div h (l) =0 in R 3 , 

Vp (l) + C h {l) = in R 3 , 



h {i) = (a on E , 



and 

dbfH (l) = 



in R 3 



VP W +CH {1) =0 in 



»3 



if (l) = e; x x on E . 

We will show that the combinations 



^fehW+WiJlW], P = J2foP (i) +"i p(i) ]+9-x, (34) 



for a suitable choice of the vectors £ and lj, solve the steady free fall problem. 
Theorem 4.2. The differential problem (27) -(32) admits a solution (u,p,£,tv,g). 
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Proof. It is straightforward to check that the fields u and p defined by (34) satisfy equations (27), 
(28), and (29). Equation (32) implies that u> = Xg for some A 6 R, and equations (30) and (31) 
reduce to the following algebraic system in the six scalar unknowns A, and g: 

K£ + ASo = m e g , 

(35) 

C£ + XBg = -m c (r X g) , 
where the matrices K, S, C, and B arc defined by 



Kji := - lim / l{h {l \p {i) )n ■ ej , 

Jav r (T.) 

Sji :=-Hm / T(H®,pV)) n .ei, 
Cji := — lim / xx T(h^' ,p^)n ■ e, , 

JdV r (T.) 

Bji := - lim / xx T(i? (l) , P {i) )n ■ ej . 



ldV r {T.) 

It is now clear that the steady free fall problem admits a solution if and only if the algebraic 
system (35) admits a solution, and the latter fact is related to the properties of the matrix 

'K S N 
C B y 

By the linearity of (33), the instantaneous energy dissipation rate E of the flow which takes the 
value £ + us x x on E is given by (see also [10]) 



E = 



c b) it) ' (t 



Now, the constitutive prescription for the Cauchy stress tensor has been chosen to be thermo- 
dynamically consistent, thus implying that E be positive for any non- vanishing flow; hence, the 
matrices A, K, and B turn out to be positive definite and, in particular, invcrtiblc. 
We can now transform (35) into 

£ = K-Ve9 + AS 9 ), (36) 

Fg := (CK^S - B)- 1 (?7i e CK- 1 9 + m c r x g) = Xg . (37) 

Since A is non-singular, the matrix F is a well-defined 3x3 real matrix and it has at least one real 
eigenvalue. Such an eigenvalue A, the associated unit eigenvector g and £ calculated as in (36), 
together with the fields u and p introduced in (34), furnish a solution for equations (27)-(32). □ 

4.2. Nonlinear free fall. To solve the differential problem (21)-(26) we employ a weak (varia- 
tional) formulation. We need to introduce the space X of virtual velocities (test fields). We first 
consider the space 

V := {v e C^(IR 3 ;R 3 ) : divt; = and v(x) = q + 1/> x x in V r (E), 3 r > 0} , 

where r/, i\) <G M 3 represent the virtual vectors associated with £ and u, respectively. Now, denoting 
by || • || 2 the norm in L 2 (R 3 ;M 3 ), the space X is the completion of V with respect to the norm 

|2 <-, no vr II 2 , „2 II a 1 1 2 



V 



x .-2\\SymVv\\ z 2 +e 2 \\Av\ ]2 



It is easy to show that, for any v G X, we have J R3 v = 0. Indeed, take any a £ M 3 and any 
v £ V: in particular v has compact support, and there exists a ball Br of radius R large enough 
to ensure that 





/ a ■ v = 


L 


7k 3 


Jb r 





The result can be extended to all of X by a density argument. 



8 



G. G. GIUSTERI, A. MARZOCCHI AND A. MUSESTI 



We multiply (22) by v 6 X, and integrate over R 3 \V^.(E) with r > sufficiently small, obtaining 



Re 



/ {[(u — U) ■ V]tt + u> x u} ■ v = / divT (u,p) ■ v + / g v 

jR 3 \V r (S) JR 3 \Vr(S) ilR3\K(E) 

= -2 / Sym Vu-Vv-i 2 Au ■ Av 

- T(u,p)n-v-£ 2 Au®n-Vv + g-j v. (38) 

ave 

/ Au®n-\7v=l div(Vu T Ait) = I VAu-Vv, 

JaVrCs) Jv r cs) Jv r (T.) 



and, letting i — > 0, 



lim / T(w,p)n • i> = lim T7 • / T(tt,p)n + if? ■ x x T(u,p)n 

r ^°Jay r (s) r ^°\ JdVrCs.) Jav r (Ti) J 

= T] f + ^ t. 

By construction Vu € i/ 1 (M 3 ; Mata(R)) and, being u a solution belonging to H 2 (M. 3 ; R 3 ), we have 
VA« e i? _1 (M 3 ; Mat 3 (R)), so that VAm • e L 1 (R 3 ;R). The latter property guarantees that, 
by taking the limit r — > in (38), we obtain 

Re / {[(u - U) ■ V]w + wx«}-d = 

Jr 3 

= -2 1 Sym Vm • Vv - I 2 I Au ■ Av - rj ■ f - ip • t ; 

JR 3 JR 3 

by substituting equations (24)-(25) into the preceding expression, we can write 

Re / {[(u — U)- V]w + uj x u} ■ v + Re(mu: x £) • 77 + Re[u) x (Ju;)] ■ V = 

Jr 3 

= -2/ Sym Vu • - £ 2 / Au ■ Av + m e g ■ rj - m c (r X g) ■ ip , (39) 

for any which represents the variational form of equations (21)-(25). 

While proving our main result we will make use of the following lemmas. The third one is 
proven in [17, Lemme 4.4]. 

Lemma 4.3. There exists a constant K > 0, depending only on £, such that, for any v G V , 

\ V \ + \ V + iPxr\ <K\\v\\ x , (40) 
where rj + xp x x is the value taken by v in a neighborhood ofT,. 

Proof. Let us first observe that, by Korn's inequality and integration by parts, \\-\\ x is equivalent 
to the norm induced on X by the standard norm in i? 2 (R 3 ;R 3 ), and that standard embedding 
theorems provide a constant K\ > such that WvW^ < K\ \\v\\ x for any v <E V. We then have 



H 1 (S)|t 7 + i/) x r| 



/ f] + ip x xdV}{x) 



and it is now enough to show that \q\ < WvW^ to prove the assertion with K = 2K±. 

By noting that the velocity of the center of mass belongs to the convex hull in R 3 of the set of 
velocities on E, we conclude that 

|f?| < sup \rj + iP x x\ < WvW^ , (41) 

a)ES 

and the proof is complete. □ 
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Lemma 4.4. There exists {<^>j}jgN C V whose linear hull is dense in X , and such that, for any 
i,k£N, 



2 / Sym V0 4 • Vcb k + t / A0,, • A0 fc = 5 lk , 

JR 3 JR 3 

and, given v £ V and e > 0, there exist m £ N and bi,...,b m el smc/i i/iai 

m 



i = l 



< e . 



(42) 



C 2 (R 3 ;R 3 ) 



Proof. Let be the completion of V in the norm 



|a|=3 

where |a| denotes the length of the multindex a. Clearly, If is a separable subspace of X, and V 
is a separable subspace of W, hence there exists a {y JieN C V dense in W. Since we are dealing 
with separable Hilbert spaces, we can obtain from {<Pj}ieN a set {</>j}ieN which is dense in X and 
orthonormal with respect to the scalar product on X which induces \\-\\ x - The approximation 
property (42) descends now from Korn's inequality and standard embedding theorems. □ 

Lemma 4.5. Fix R > and let Br denote the ball in K" of radius R and center 0. Let 

T : x S 2 -¥ W , 
t : E7r~ x S 2 -> M 3 

&e continuous maps. If 

V(c,g)£dB R x S 2 :T(c,g)-c>0, and V (c, g) £ B R x S 2 : r(c, g) • g = , 
t/ien i/iere exists (c, g) £ -Br x S" 2 smc/i i/iai 

T(c, g) = and r(c, g)=0. 
Theorem 4.6. TTie differential problem (21)-(26) admits a solution (u,p,£,u),g). 
Proof. We consider, for any n £ N, the truncated expansions on the basis furnished by Lemma 4.4 



4=1 



where <p i = r] i + ip i x a; on S, and apply Galerkin's approximation to equation (39). 
For any n £ N, we consider the maps defined by 



T(c,g) := g x w r 



and 



[T(c,g)] k ■■= Re / {[(u n - £„ - w„ x x) • V]«„ + u>„ x m„} • d> fc 

+ Re(mu; n x £„) • rj*. + i?e[u;„ x (J«„)] • Vfe 
+ 2 / Sym Vm„ • V0 fe + ^ 2 / Au n • A0 fe - m e g • r] k + m c (r x g) • , 

JR 3 JR 3 

where we denoted by c the rt-dimensional vector of the coefficients of the truncated expansion 
which defines u n . 



10 



G. G. GIUSTERI, A. MARZOCCHI AND A. MUSESTI 



It is immediate to verify that r and T are continuous maps and that, for any (c,g) S Br x S 2 
we have t(c, g) ■ g — 0. Moreover, 

T(c,g) ■ c = Re / {[(u n - £ n - w„ x x) ■ V]tt„ + w n x «„} • u n 

+ Re(mu: n x £J • £„ + fle[w„ x (Jw n )] • w„ 
+ 2/ Sym Vit„ ■ V«„ + / Aw„ • Aw„ - m e g • + m c (r x gr) • w„ 

= 2 / I Sym Vm„| 2 + I 2 / |Am„| 2 - m e g • + m c (r x g) • «„ . 

JR 3 JR 3 

In the previous equality we used the fact that 

/ [(«n - £„ - 0J n X x) ■ V]tt n ' «„ = , 

JR 3 

which is easily proved, since divu n = 0, via tensorial identities and integration by parts. Now, by 
using Lemma 4.3, there exists a constant C > 0, independent of n, such that 

T(c,g) - c = - m e g ■ £ n + m c (r x g) • u n 

= R 2 - (m e + m c )g ■ £„ + m c g ■ (£„ + u n x r) 

> i? 2 - ro(|£ n | + |£ fl + w„ x r|) > i? 2 - Ci? > , (43) 

for i? large enough and whenever (c, g) € c>-B,r x S 2 . Hence, Lemma 4.5 guarantees the existence, 
for any n € N, of a pair (c,g n ) that produces a solution (u n ,^ I1 ,u„,g n ) for the approximate 
problem given by 

g n x u>„ = , 

and 

iZe / {[(«„ - £„ - u„ x x) • V]u„ + u>„ x m„} • </> fc 

JR 3 

+ Re(mu: n x £J • ?7 fe + i?e[w„ x (Jw n )] • Vfe 
+ 2 / Sym Vw„ • V0 fc + I 2 / A« n • A0 fe - ?n e g„ • r; fc + m c (r x t/„) • ip k = , 

JR 3 JR3 

for any k £ N. 

By taking into account (43), we see that the sequence {u n } raS N is bounded in X, and hence 
weakly convergent to some u £ X; moreover, u n — > u strongly in L 2 (R 3 ; M 3 ), and, clearly, £ n — > £, 
uj n — > u>, and g n — > g in M 3 . Since the nonlinear term 

is a compact operator from X into its topological dual (see [8, Theorem 1]), the limits of the 
previous sequences are such that 

g x w = 0, 

and 

i?e / {[(w — £ — w x x) ■ V]u + lj x u} ■ (fi k 

Jr 3 

+ Re{mui x £) • rj k + Re[uj x (Ju>)] • xp k 

+ 2 Sym Vm • + f 2 / Au ■ A<fi k - m e g ■ r/ k + m c (r x g) ■ xp k = , 

for any fc s N. Finally, by suitable linear combinations of the fc 's and the approximation property 
granted by Lemma 4.4, we can show that (u,£,u,g) is a solution of equations (26) and (39) for 
any v G X. The pressure field p needed for the solution of the problem (21)-(26) can be recovered 
as the Lagrange multiplier of the divergence- free constraint imposed on the space X. □ 
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5. Unsteady free fall 

We now consider the time-dependent differential problem given by (9)-(16). The existence of 
a weak solution can be easily proved by standard methods, as in the classical situation considered 
in [4], but, in the present context, the additional regularity provided by the hyperviscous term 
allows to prove also uniqueness and continuous dependence on the initial data of the solution for 
the free fall problem. 

Theorem 5.1. Let H be the completion of V with respect to ||-|| 2 . Given any Uq £ H and any 
g G S 2 , there exist a unique u belonging to the space 

X := L 2 ([0,T};X)nC a ([0,T};H)r\H 1 ([0,T]]X') 

with u{x,t) = £(t) +u)(t) xx onEx [0, T] and u(x, 0) = u (x), and a unique G G 7J 1 ([0, T]; 5* 2 ) 
with G(0) = g such that 



du 
~dt 



Re[(u - U) ■ V]n + Reu) x u 



I— 

dt 



- 2 



JR 3 



Re(muj x 



Sym Vm • Vv - t 



I duj 

~dt 

T 



Re(u x (Ju>)) 
Am • Av 



f (m e G ■ r] - m c (r X G) ■ ip) , (44) 
Jo 



for any and 



dG 

dt 



= Re(G x u) 



(45) 



Proof. The existence of a weak solution to the free fall problem, as defined by (44)-(45), can be 
obtained by standard techniques exploited in the theory of the classical Navier-Stokes equation 
(see e.g. [15]). We then proceed to prove the uniqueness of the solution. We take two solutions, 
represented by (u\, £ 1; Wi, Gi) and (U2, £ 2 i ^2, G2), and denote their difference by (u, £, w, G). 
Since u G X, both solutions satisfy (44) with v = ii and (45) tested on G; by considering the 
difference of such equations and repeatedly applying Young's inequality one easily obtains, for any 

te[0,T], 



\ m \m 2 + \^{t)-"{t) + \\G{t)\ 2 



u 



X 



Re 



= -Ref f ([(m-L0-V]w 2 + (wxm 2 )}-m 
J |(mw x ^ 2 ) • ^ + (uj 2 x Jo>) • d) - (G 2 x ib) ■ g| 



< 



u 



x 



Ci 



m e G ■ I — m c (r x G) ■ Cj 



\u 2 \\l\\u\\t + C 2 



m 2 + c, / ig 



Denoting by ji the least eigenvalue of the inertia matrix J, which is positive definite, we finally 
have 



\u(t)\\ 



|^| 2 +J ^(0| 2 + |G(i)| 2 <C< 5 



u 



m\£\ 2 + n \u>\ 2 



and, by Gronwall's lemma, we can conclude that (ii,^,a>,G) = (0,0,0,0) for almost every t G 
[0, T], thereby showing that the solution is unique. □ 
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